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Abstract 

The ROSAT all-sky survey discovered several 'shrapnels', showing boomerang 
structures outside the Vela supernova remnant. We observed shrapnel D with the 
XMM-Newton satellite. There is an X-ray bright ridge structure in our FOV running 
from north to south. Applying the VNEI model to X-ray spectra of various regions, 
we find that the plasma in the eastern part from the X-ray ridge is significantly dif- 
ferent from that in the western part. The X-ray spectra in the western part can 
be represented by a single-temperature component. The abundances of heavy ele- 
ments are almost uniform, whereas they are heavily overabundant, except for Fe; the 
relative abundances to the solar values are O ~ 5, Ne ~ 10, Mg ~ 10, Fe ~ 1. This in- 
dicates that shrapnel D originated from the ejecta of the supernova. We find that the 
plasma in the eastern part from the ridge consists of two components with different 
temperatures; the hot component comes from the ejecta, while the cold component 
comes from the interstellar matter. These two components are considered to be in 
contact with each other, forming a contact discontinuity. Around the northern part 
of the contact discontinuity, we find wave-like structures of which the typical scale 
are comparable with that of the Rayleigh-Taylor instability. 

Key words: ISM: abundances — ISM: individual (Vela Supernova Remnant) - 
supernova remnants — X-rays: ISM 



1. Introduction 

The Vela supernova remnant (SNR) is a middle-aged SNR that exploded about 10000 
years ago (Taylor et al. 1993). It clearly shows a circular structure with a diameter of 8°. 3 
(Aschenbach et al. 1995). At a distance of 250 ± 30 pc (Cha et al. 1999), this corresponds 
to 36 pc. The Vela SNR is so close to us that it can be an ideal candidate for resolving fine 
structure in X-rays. Aschenbach et al. (1995) also discovered 'shrapnels', boomerang structures 
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from 'A' to 'F', outside the main shell. The opening angles of the shrapnels suggest supersonic 
motion in tenuous matter. The center line of each boomerang structure can be traced back 
to the center of the SNR, which is close to the Vela pulsar. Therefore, they are considered to 
originate from fragments of supernova (SN) ejecta that are beyond the position of the main 
blast wave. There is an alternative explanation that the features are "break-outs" of shock in 
which inhomogeneities in the ambient medium cause the shock to be non-spherical. 

Vela shrapnel A was observed with ASCA (Tsunemi et al. 1999). It was clarified that 
the abundance of Si is about 10-times higher than that of O. Therefore, they concluded that 
shrapnel A was an explosion ejecta from a Si-rich layer of a progenitor star. If it is ejecta 
of a SN explosion, the interstellar matter would be swept up in the leading edge while the 
ejecta material would be peeled off in the trailing edge. However, the spatial distribution of Si 
could not be measured due to the poor spatial-resolving power of ASCA. In order to investigate 
the spatial structure of shrapnel A, Miyata et al. (2001) observed it with the Chandra X-ray 
Observatory. The Chandra image reveals a bright X-ray region at the head position of shrapnel 
A and a fainter extended tail. They confirmed an overabundance of Si in the head region, but 
could not obtain a strong constraint on the parameters of the tail region. 

Shrapnel D, the eastern limb of the Vela SNR, is the closest to the main shell and the 
brightest in X-rays of the six shrapnels. The optical nebula RCW 37 (Rodgers et al. 1960) 
lies along the outer edge of shrapnel D, which clearly shows that shrapnel D is now interacting 
with an interstellar cloud. Sankrit et al. (2003) observed this optical filament and found that 
shrapnel D is a bow shock propagating into an interstellar cloud with normal abundances. In 
order to understand the nature of shrapnel D, Plucinsky et al. (2002) observed it two times 
with Chandra, one at the head of shrapnel and the other in the trailing "wake" . If the X-ray 
emission associated with the shrapnel D is produced by shock-heating of the ambient medium 
by supersonic motion of the ejecta, one would expect the abundances to be enhanced at the 
head region and close to the interstellar medium (ISM) value at the wake region. They found 
that the spectra from the different locations in the shrapnel are remarkably similar to each 
other. They did not show the abundance variations that might be expected from a fragment 
of ejecta. However, they could not obtain strong constraints on the elemental abundances: 
the O and Ne abundances vary from 0.5 to 5.0 x solar and 1.6 to 6.4 x solar, respectively. 
They therefore concluded that the origin of shrapnel D is more consistent with a shock breakout 
hypothesis than with an ejecta hypothesis, although they cannot rule out the ejecta hypothesis. 

What is indeed the origin of shrapnel D? If it is a shrapnel of ejecta, how is the dis- 
tribution of the element abundances in it? In order to answer these questions, we observed 
shrapnel D with the XMM-Newton satellite. We here present the results from the observation. 
Then, we performed a spatially resolved spectral analysis, and showed the results of an in-situ 
observation of interaction zones between shrapnel D and the interstellar cloud. 
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2. Observations 



We observed shrapnel D on 2002 May 2 with the EPIC (European Photon Imaging 
Camera) instruments on board the XMM-Newton satellite. The observation ID is 0136010101 
(PI H. Tsunemi). Since the shrapnel is larger than the field of view (FOV) of the EPIC (30'), we 
selected the brightest region in the shrapnel. Figure 1 shows a ROSAT PSPC image of shrapnel 
D with an overlay of the FOV of our observation. The EPIC MOS cameras were operated in 
the standard full-frame mode and the EPIC PN camera in the extended full-frame mode. Each 
camera was operated using a medium filter. Since we screened the data by rejecting the high- 
background periods, the effective exposure time was reduced to ~ 21 ks for MOS-1 and MOS-2 
and to ~ 14 ks for PN. We selected X-ray events corresponding to patterns 0-12 for MOS and 
patterns 0-4 for PN, respectively. 

Since shrapnel D is larger than the FOV, we could not subtract the background events 
from the same data. Therefore, we used an observation (ID 0147511701) of the Lockman hole 
as background data. 

3. Image Analysis 

Figure 2 shows an exposure-corrected MOS-1 + MOS-2 image of shrapnel D in the energy 
range of 0.3-2 keV with an overlaid optical contour map. We extracted the exposure map using 
SAS v5.4.1(the EEXPMAP command), and smoothed the image by a Gaussian of a = 6". We 
found a bright X-ray ridge structure running from north to south in the region east of the FOV. 
The X-ray emission from the eastern part of the ridge abruptly weakens to the background level, 
while that from the western part of the X-ray ridge gradually weakens. There is also an optical 
bright ridge structure, a pencil nebula RCW 37, running parallel to the X-ray ridge at about 3' 
in the east. This suggests that shrapnel D is now interacting with an interstellar cloud (Sankrit 
et al. 2003). 

4. Spectral Analysis 

We investigated spectral variations from various regions in detail. According to the X- 
ray ridge structure, we divided our FOV into many small rectangles (2' x 4', a side with 2' is 
perpendicular to the ridge and that of 4' is parallel to that). The scale of the rectangle is larger 
than the half-power beam width of XMM-Newton. The scale was selected such that we could 
obtain sufficient photon statistics in each rectangle. For the purpose of more detailed investi- 
gations, we arranged rectangles such that each rectangle overlapped each other by half of its 
size, obtaining 243 spectra in total. Figure 3 shows those regions overlaid on a MOS-1 + MOS-2 
image that is the same as figure 2. In the spectral analysis, we included the uncertainties among 
CCDs by 5% on the model in quadrature (see, e.g., Nevalainen et al. 2003; Kirsch 2004). We 
found that the spectra from the western part of the X-ray ridge were quite different from those 
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from the eastern part (the boundary is indicated in figure 3). Therefore, we divided the data 
into two regions by the X-ray ridge. We present our results in the following subsections. 

4-1- Western Part of the X-Ray Ridge 

Figure 4 shows an example spectrum from the white rectangle in figure 3. We can 
see prominent O VII triplets at ~ 0.57keV, O VIII Lya at ~ 0.65 keV, Ne IX triplets at 
~ 0.91 keV, Ne X Lya at ~ 1.02 keV, Mg XI triplets at ~ 1.35 keV lines. These features 
are clearly seen in all of the spectra from the western part of the X-ray ridge. We applied 
an absorbed non equiliblium ionization model (the VNEI model in XSPEC v 11.3.1) to each 
spectrum (Hamilton et al. 1983; Borkowski et al. 1994, 2001; Liedahl et al. 1995). The free 
parameters in our analysis are the electron temperature, kT e ; the ionization parameter, r; the 
emission measure (hereafter EM; EM = / n e riidl, where n e and n\ are the number densities of 
electrons and ions, respectively, and dl is the plasma depth); the abundances and the column 
density, Ah, of the absorbing foreground material. Above, r is the electron density times the 
elapsed time after shock heating. Only the abundances of C, N, O, Ne, Mg, Fe are free, and 
the others are frozen to the solar values (Anders, Grevesse 1989). We set the abundances of 
C and N equal to that of O. All of the spectra are represented by a single-temperature VNEI 
model. 

The reduced \ 2 is 1—1-5 for about 85% of the region to the west of the X-ray ridge. 
The remaining (only 15% in the west of the ridge) regions, which are located near the ridge 
and have a high surface brightness, show about 1.5-1.8. These values are far from acceptable 
from a statistical point of view. These regions are very close to the interaction region with 
the interstellar cloud. We could not obtain better fits by adding extra components there. 
Therefore, our model is too simple to reproduce the entire region. A more advanced model will 
be needed to obtain a better fit, particularly near the X-ray ridge. In table 1, we show the 
best-fit parameters of the spectrum from the white rectangle. 

Figure 5 shows the values of kT e and a logr contour map in the western part of the 
X-ray ridge overlaid on the MOS-1 + MOS-2 image. The value of kT e shows the highest 
(0.32 ±0.01 keV) at the ridge and a gradual decrease in the west down to 0.25 ±0.02 keV. The 
spatial variation of logr shows an anti-correlation with that of kT e . Most of the regions, having 
relatively low kT e , show r > 10 12 scm -3 , which indicates a collisional ionization equilibrium 
condition. However, the regions near the X-ray ridge, having a relatively high kT e , show that r 
is (4.8-7.2) xlO n scm -3 , which indicates a nonequilibrium ionization condition. Figure 6 shows 
contour maps of metal abundances. They are heavily overabundant, except for Fe, at any re- 
gions in the map; the values of the abundances relative to the solar values are O ~ 5, Ne ~ 10, 
Mg ~ 10, Fe~ 1. These values support the idea that shrapnel D originates from the explosion 
eject a. 
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4-2. Eastern Part of the X-Ray Ridge 

A weak X-ray emission comes from the eastern part of the X-ray ridge. We found that 
the spectra from the eastern part of the X-ray ridge were quite different from those of the 
western part. We divided the eastern part into two regions. 

4-2.1. Southeastern region near the optical filament 

We found that the spectrum from the region near the optical filament could be well fitted 
with a single component with solar abundances. It shows a relatively low kT c (~0.15keV). 
Sankrit et al. (2003) found that optical emission comes from a shock-excited cloud with solar 
abundances. Therefore, we conclude that the \ow-kT c component represents plasma originating 
from the cloud. 

We also found that the spectrum varies between the optical filament and the X-ray 
ridge. We analyzed the spatial variation of the spectra between them. We extracted spectra 
from rectangular regions (1' x 8') drawn by the white line in figure 7 a, each of which overlapped 
adjacent regions by the half width. In this way we extracted 5 spectra. We found that the 
spectra were not represented by a single-temperature VNEI model. Therefore, we added an 
extra component with different kT c . Due to the insufficient photon statistics of the low- kT c 
component above the O-line energy, we assumed solar abundances for all elements, except O. 
The extra component improves the fit (e.g., reduced \ 2 is improved from 1.7 to 1.2), which is 
in contrast to that of the western part. The kT e values of these two components are 0.15 keV 
and 0.35 keV, while the intensity fraction varies from east to west. The abundances of the high 
kT e component are similar to those of the western part of the X-ray ridge, which shows that 
it comes from the ejecta. The abundance of O and kT e for the \ow-kT e component are similar 
to those for the cloud component. Therefore, the low-A;T c component is considered to originate 
from the cloud. The fraction of the cloud component decreases toward the X-ray ridge. We 
show the variation of EM corrected by the interstellar absorption as a function of the location 
in figure 8. The EM of the cloud component decreases toward the X-ray ridge, while that of 
the ejecta increases toward the same direction. It is reasonable to consider that this feature is 
due to the projection effect, as looking obliquely between the cloud and the ejecta. We consider 
that the boundary region between these two plasmas is a contact discontinuity. 

4-2.2. Northeastern region 

In the X-ray surface brightness map, we found wavelike structures at the northeastern 
region of the FOV in figure 7 a. Figure 7 b shows a band ratio map of the O VIII band (624- 
684eV) / the O VII band (545-603eV), in which the wavelike structures are prominent. The 
arrow in figure 7 b shows A D , a typical length of the wavelike structure, to be ~ 0.3 pc. 

In order to investigate the feature of the wavelike structure, we extracted spectra from 
regions indicated in figure 7 b. The circular regions with the thick solid line and dashed line 
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are responsible for the regions where the O VIII band is dominant, while those with thin line 
are responsible for the regions where the OVII band is dominant. All of the spectra are 
represented by a two-component VNEI model, just as those applied in the southeastern region. 
It is reasonable to consider that the low-temperature component represents an interstellar cloud, 
while the high-temperature component represents the ejecta. Then, the cloud dominates in the 
circular regions with thin line and the ejecta dominate in the circular regions with thick line. 
We consider that there is a contact discontinuity in those regions. Figure 9 shows the spectrum 
from the dashed circle indicated in figure 7 b and the best-fit parameters are given in table 1. 
The value of the plasma depth of the ejecta comes from the assumption that the circular shape 
in figure 7 is due to the spherical shape. This result suggests that the surface of the contact 
discontinuity has wavelike structures with a typical length, A Q . 

5. Discussion 

X-ray spectra in the western region of the X-ray ridge can be represented with a single- 
temperature VNEI model. The abundances of heavy elements, except for Fe, are significantly 
higher than those of solar values in any region in the western part of the X-ray ridge: O ~ 5, Ne 
~ 10, Mg ~ 10, Fe ~ 1. The absolute abundance obtained by the spectral fits depends on the 
model employed, whereas the relative abundances among the heavy elements are robust. The 
relative abundances, (O/Fe, Ne/Fe and Mg/Fe) are ~ 5, ~ 10 and ~ 10, respectively. High-Z 
elements, like Ar, Ca, and Fe, are generated dominantly in a Type-la SN (Nomoto et al. 1984), 
while low-Z elements, like O, Ne, Mg, and Si, are generated in a Type-II SN (Thielemann et al. 
1996). The relative abundances indicate a Type-II SN origin rather than Type-la. Since Vela 
SNR is believed to be a Type-II SN, we conclude that shrapnel D is ejecta originating from an 
explosion of Vela SN. 

We estimated the total mass of the ejecta of Vela shrapnel D. We assumed that shrapnel 
D has a conical structure with a diameter 80' and a height of 45'. Using a distance of 250 pc, 
and the assumption of uniform electron density of ~ 0.2 cm -3 , we can estimate the total mass 
to be about 1O _1 M . 

We found that shrapnel D is ejecta coming from O, Ne and Mg-rich layer of a progenitor 
star. This is contrast to the fact that shrapnel A comes from the Si rich layer (Tsunemi et al. 
1999, Miyata et al. 2001). These findings suggest that shrapnel D comes from an outer layer of 
the progenitor star than shrapnel A. We compare the present distance from the center of Vela 
SNR to shrapnel D with that to shrapnel A. Aschenbach et al. (1995) inferred the temperature 
of shrapnel D to be 0.34±0.07keV from the opening angle of the feature and the temperature 
of the ambient medium. The opening angle was estimated based on the assumption that the 
shrapnel is moving in the plane of the sky. The measured temperature (0.32±0.01 keV) is 
in good agreement with that inferred from the geometry, which supports that shrapnel D is 
moving in the plane of the sky. Therefore, the actual distance of shrapnel D from the center 
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of Vela SNR is shorter than that of shrapnel A. Since shrapnel D is considered to have come 
from the outer layer of a progenitor star than shrapnel A, it suggests that a part of the inner 
layer has a higher initial velocity than that of the outer layer. 

In the southeastern region of the FOV, there is an optical filament, a pencil nebula, RCW 
37, running parallel to the X-ray ridge about 3' to the east. We confirmed that the emission 
from the optical filament region comes from a shock-excited cloud (Sankrit et al. 2003). We 
found the spectrum variation in the region between the X-ray ridge and the pencil nebula. 
An extra component is needed to fit the spectrum. The low-temperature component shows 
solar abundances, while the high-temperature component shows high metal abundances. The 
fraction of the low-temperature component gradually decreases toward the X-ray ridge, while 
that of the high-temperature component increases. Therefore, the low-temperature component 
is considered to come from the cloud and the high-temperature component is considered to 
come from the ejecta. 

We found wavelike structures in the northeastern region of the FOV in the X-ray surface 
brightness map. This structure is more apparent in the band ratio map between the OVIII 
band and the O VII band. Furthermore, our spectral analysis indicates that a contact dis- 
continuity between the ejecta and the interstellar cloud has wavelike structures. The wavelike 
structures remind us of the Rayleigh-Taylor (R-T) instability, which occurs where a heavy 
fluid is accelerated by a light one (i.e., the pressure of the light fluid is higher than that of the 
heavy one). In our case, the temperature of the measured ejecta (T ej -) exceeds that of the cloud 
(T c i). Therefore, if this is the case, the heavy fluid must be the cloud and the light one must 
be the ejecta. The typical scale of R-T instability is given by (Velazquez et al. 1998) 

5 \ 1/3 



A = 0.39 x ( Tj 9 ) pc, 
\g 2 an{) 



where 




^ 2 10 2 cms 2 ' \n c \ + n Cj 
Here, g e g is the effective gravity (inertial force on a frame of the contact discontinuity), 
7i c i, the number density of the cloud and n e j that of the ejecta. We assumed the cloud depth 
to be 0.73 pc (= 10' ) in order to estimate n c \. The value of g c s is given by (1/ p m )(dP/dr) 
(Ebisuzaki et al. 1989), where dP/di is the pressure gradient and p m is the mean density 
between the cloud and the ejecta. We used the approximation that dP/di ~ AP/X, where AP 
represents the pressure difference between two fluids. We consider that the temperatures of the 
cloud and the ejecta come from shock heating. After shock heating, the ion temperature would 
not be equal to that of the electron. Then, the electron and ion temperatures vary smoothly to 
reach the thermal equilibrium. The postshock electron and ion temperatures equilibrate rapidly 
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after a slow shock, even when the shock is collisionless (Ghavamian et al. 2001). Therefore, we 
assumed thermal equilibrium between the ion and the electron, P\ = P c (Pi,P e are the ion and 
the electron pressure, respectively), because the shock speed is only about 150 kms" 1 (Sankrit 
et al. 2003) in our case. Based on the assumption of P\ = P c , the pressure difference between 
the cloud and the ejecta becomes 2 x AP e (AP e is the electron pressure difference between the 
cloud and the ejecta). Using the values obtained by the spectral fitting (from the dashed circle 
shown in table 1), we obtained A as 



where 



0.186/ V3.7/ 
P 



a ~ 0.2 x ( 5/2 r ^ 1/2 ( % v 1/2 f^r /2 f^v 1 pc 



"ml 



1.67 x 10- 24 gcm" 3 ' 

Pro = AP 

1 x 10 _10 ergcm~ 3 

We found that A is close to A Q . Since the value of n strongly depends on the plasma 
depth, the value of A also depends on the plasma depth. It is highly probable that the wavelike 
structures are caused by the R-T instability. 

In the southeastern region of the FOV (around the optical filament), we cannot see 
wavelike structures. In this region, we find n e j < n c \, and the pressure of the ejecta is lower 
than that of the cloud, where the contact discontinuity is stable. 

6. Conclusion 

We observed shrapnel D with the XMM-Newton satellite. We found an X-ray bright 
ridge structure running from north to south. We found that the spectra from the western part 
of the X-ray ridge could be represented by a single kT c component with high metal abundances: 
O ~ 5, Ne ~ 10, Mg ~ 10, Fe ~ 1. This clearly shows that the origin of shrapnel D is ejecta 
from Vela SN. The spectra from the eastern part of the X-ray ridge were represented by two kT e 
components. The low-temperature component having solar abundance is considered to come 
from the interstellar cloud, while the high-temperature component is considered to come from 
the ejecta. This is the place where the ejecta are interacting with the interstellar cloud. In the 
southeastern region of the FOV, it is reasonable to think that we see the interaction obliquely 
due to the projection effect. In the northeastern region of the FOV, it is highly probable that 
we see the wavelike structure of the R-T instability. 
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Table 1. Spectral-fit parameters.* 



Parameter White rectangle in figure 3 Dashed circle in figure 7 

Ah KPnn 2J+^ mT 7 ^ 

Low temperature component 

kT, [keV] 0.16±g;g| 

o(=c=n) o.ii 

T[scm- 3 ] >3x 10 11 

EM[cm- 5 ]t (13±2)xl0 17 

Plasma depth [pc] 0.73 

Electron density [cm~ 3 ] 

Electron pressure [10~ 10 crgcnT 3 ] 1.99±0.22 

High temperature component 

fcT [keV] 0.30±0.01 0.38±0.01 

0(=C=N) 6.4±^ 4.5±°;? 

Ne 13.7±°; 8 3 11.51° j 

Mg 14±2 8.9^2:? 

Fe 1.41°;? 1.00t°;i* 

r[scm" 3 ] (6.1 ±0.2) x 10 11 (1.9 ± 0.06) x 10 11 

EM[cm" 5 ]t (2.60 ± 0.02) x 10 17 2.66±g;g| x 10 17 

Plasma depth [pc] 2.2 0.2 

Electron density [cm" 3 ] 0.197±0.001 0.627±0.007 

Electron pressure [10- 10 crgcnr 3 ] 0.96±0.06 3.86±0.25 

X 2 /d.oi 436/303 362/257 

* Other elements are fixed to those of solar values. The values of abundances are multiples of solar value. 
The errors are in the range A\ 2 < 2.7 on one parameter. In the fits, 5% systematic errors are included. 
tEM denotes the emission measure J n e nidl. 
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figure 1 . eps 



Fig. 1. ROSAT PSPC logarithmically scaled image of the Eastern Limb of the Vela SNR. The black circle 
indicates the FOV of our XMM-Newton observation. 



f igure2 . eps 



Fig. 2. EPIC MOS-l+MOS-2 logarithmically scaled image in the energy range 0.3-2.0keV. The data 
have been smoothed with a Gaussian of a = 6". The contours show a linearly scaled optical intensity map. 
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f igure3 . eps 



Fig. 3. Same as figure 2 with an overlaid all the regions where we extracted spectra. The spectrum from 
the white rectangle is shown in figure 4. 




0.5 1 Z 

Energy [keV] 

Fig. 4. X-ray spectra extracted from the white rectangle in figure 3. The best-fit curves are shown with 
solid lines and the lower panels show the residuals. 
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f igure5_a_r200 . eps 



f igure5_b_r200 . eps 



Fig. 5. Linearly scaled contour map of kT e and logr in the western part of the X-ray ridge overlaid on 
an X-ray image the same as figure 2. The values of kT e are in units of keV. The data were smoothed with 
a Gaussian of a = 4'. 



f igure6_a_r200 . eps 



f igure6_b_r200 . eps 



f igure6_c_r200 . eps 



f igure6_d_r200 . eps 



Fig. 6. Linearly scaled contour map of each metal abundance overlaid on an X-ray image the same as 
figure 2. The values shown in each figure are those relative to the solar abundance. The data were 
smoothed with a Gaussian of a = 4'. 
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f igure7_a050710 . eps 



f igure7_b050710 . eps 



Fig. 7. (a) Same as figure 2. The white rectangles adjacently lined up are the regions where we extracted 
the spectra. The large upper-left region indicates the region where we extracted the band ratio map. (b) 
EPIC (MOS-1 + MOS-2 + PN) band ratio map (OVIII band / OVII band). The circles are the regions 
where we extracted spectra. 
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Fig. 8. Variation of EM of two components in the region between the optical filament and the X-ray ridge 
indicated in figure 7 by white rectangles. 
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Fig. 9. X-ray spectra extracted from the dashed circle in figure 7 b. The solid line is the total model, 
while the broken lines represent the individual contributions. The lower panels show the residuals. 
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